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Abstract. - A statistical mechanics expression is derived for the space-dependent dielectric 
permittivity of a polar solvent near an interface. The asymptotic behaviour of this local per- 
mittivity is calculated in the low density limit, near a planar interface. The potential of mean 
force between two ions is shown to agree with the prediction of macroscopic electrostatics for 
large separations parallel to the interface. 



Introduction. - An efficient statistical description of complex fluids and biomolecular 
assemblies generally requires some degree of coarse-graining to be tractable. A good example 
is the "implicit solvent" description of colloids or biomolecules in aqueous solution where water 
is replaced by a polarizable continuum. An important, largely unsolved, question that arises is 
that of inhomogcneous screening of solute ions near interfaces. Can the reduction of the bare 
Coulomb interactions between charges be described by a spatially varying permittivity t(r), 
and can one establish a rigorous statistical mechanics foundation of such a local permittivity, 
similar to Kirkwood's famous expression 1 relating the constant permittivity e of a bulk fluid 
to dipolar fluctuations? While many phcnomcnological prescriptions have been put forward for 
the determination of local permittivities 2 , the second question is so far largely unanswered. 
In a recent paper ^ , an attempt was made to calculate the permittivity of a polar fluid near a 
planar interface by a generalization of Onsager's cavity model 111. In this letter, we reconsider 
the problem in a general statistical mechanics framework and a number of exact results are 
established to assess the validity of the concept of a local permittivity. 

Molecular model. - We consider a classical polar fluid made up of N identical molecules 
confined to a volume V , at temperature T . The surrounding media that confine the fiuid 
may be polarisable, and are treated using macroscopic electrostatics. The molecules carry a 
permanent dipole moment p and are assumed to be non polarisable. The position and the 
dipole moment of the iih molecule are denoted collectively by i = (rt^pi). The molecules 
interact via the pair potential 

2) = 2) + e(ri2-fT)(pi-Vi)(p2-V2)G(ri,r2), ri2=r2-ri, (1) 

where the first term is a short-range interaction(^) and the second term is the dipolar in- 
teraction with a cut-off for distances ri2 = \ri2\ less than the molecular diameter a [O is 

(^)The potential Vsj:{l, 2) must decay at least as |ri2|^'' at large relative distances. In addition to short-range 
repulsion, it may include the effects of higher multipole moments of the molecules. 
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the Heaviside function]. The Green's function G{ri,r2) is the solution of the equation 
AG(ri,r2) — 47r(5(r2) for ri,r2 G V, with the appropriate boundary conditions dictated 
by the surrounding media. In the case of a planar interface between a continuous dielectric 
medium of dielectric constant e' in region z < and a polar fluid in region 2: > (see Fig. 1), 
G is given by 



G(ri,r2) 



1 



1 



1 



|ri -rsl 1 + e' In -r*|' 



(2) 



where the second term corresponds to the interaction with the image charge located at 
The presence of the surrounding polarisable media gives also rise to the one-body potential 



t/(l) = i hm (pi-Vi)(pi-V2) 



G(ri,r2)- 



1 



\ri - r2\ 

which corresponds to the interaction of a dipole with its own image(s).(^) 



(3) 




Fig. 1 - A polar fluid confined to the region z > 0. 



Polarisation in the linear regime. - When an external electric field Eo{r) — — V(/)(r) is 
applied to the fluid('^), the latter acquires a non vanishing mean polarisation given by 



P(r) = J p^{r,p)pdn 



(4) 



where p^{r,p) is the one-body density. Since the coupling energy of the i*^ molecule to the 
external field is ip{i) ~ ~Pi ■ Eo{ri), the variation 5p{l) = p^{l) ~ /C0=o(l) of the one-body 
density due to the perturbation is given in linear response by 



<5p(l) = -/3^(l)p(l)+p(l) / d2Ml,2)p(2)(-/3V(2)) 



(5) 



where 2) is the pair correlation function of the unperturbed system. Using the Ornstein- 
Zernike equation, we see that 6p(l) is linked to the direct correlation function via 



6p{l)^-Pi;{l)p{l) + p{l) / d2c(l,2)<5p(2) 



(6) 



In a polar fluid, we expect the direct correlation function to decay at large distances like 
the dipolar potential. The integral in eq. JHJl depends thus in general on the shape dV of 

(^)The factor i occurs because, in a charging process where the dipole moments of the molecules are increased 
from zero to their final value, the infinitesimal energy increment dU{l) is a quadratic function of |pi|. 
{^)EQ{r) is by definition the electric field in the absence of the fluid (the subscript stands for 'no molecules'). 



V. BALLENEGGER et al: LOCAL DIELECTRIC PERMITTIVITY NEAR AN INTERFACE 



3 



the sample. Following Ramshaw 'SI, let c(l,2) = Csr(l,2) + Cdd(l,2) where Csr(l,2) is a 
short-ranged function and Cdd(l, 2) = — /3fdd(l, 2) where 

Wdd(l,2)-e(ri2-A)(pi-Vi)(p2- V2)G(ri,r2) + yPi-p2'5(r-i2), A ^ 0, (7) 

is the electrostatic interaction between two ideal dipoles. The definition of Cdd(l,2) at short 
range is of course arbitrary, and the delta function, which appears naturally in the definition 
of the dipolar potential (see 0), is included here for convenience. The cutoff 0(7-12 ^ A) 
excludes the singularity of the dipolar potential at the origin. Integrals of Vdd(l, 2) are hence 
well defined (for finite volumes V), and the limit A — > is taken at the end of the calculation. 
Using Q and lO, we obtain 

Mi) = p(i)/5pi-£^(n) + p(i) J d2cs,il,2)Spi2) (8) 

where E{ri) = -Eo(ri) — Vi / {P{r2) ■ V2)G(ri, r2) dr2 is the average electric field at ri. 
Notice that the macroscopic field [and not the Lorentz field E\^{r) = £^(r)+47rP(r)/3] appears 
in eq. ((HJ because the delta function term provides exactly the "missing" contribution from 
the small excluded volume of radius A. The shape dependence is now entirely embodied in 
the Maxwell field E{ri) since the integrand in eq. © decays rapidly. By iterating eq. (|HJl and 
substituting into eq. (0J, we find 

P(ri) =p(l)/3pi •£;(ri)+/3 j AU^ j d2 p{l)p{2)h,,{l,2)p^{p2 ■ Eir^)) (9) 

with 

/l,,(l,2) =Csr(l,2)+ /d3c,r(l,3)p(3)/^sr(3,2). (10) 



Equation relates the polarisation to the average electric field via an integral kernel. A 
local permittivity tensor e(r) can be defined when this relation is local. In that case, 

eir) - 1 

P(r) = -^^E{r) (11) 

with 

ey(ri) = (5y +4^/3 j dr2 J (pi • e,)(p2 • e,)p{l)[6il,2) + p{2)h,,{l,2)]. (12) 

This is the general statistical mechanical expression of a local dielectric permittivity e(r). For 
points ri far from the boundaries, the fiuid is homogeneous and isotropic [p(l) — p/(47r)], and 
eq. p2|) reduces to the bulk dielectric constant of the fluid: 



e- 1 1 



3y 4tt 



dr2 / driidri2Pi ■ P2 



6(1,2) + £-hr\l,2)], y='-^. (13) 



The present calculation follows the standard approach |7] to the molecular theory of the 
dielectric constant based on the decomposition of the direct correlation function into a short- 
range and a dipolar part.('*) The condition that the electric field must not vary appreciably 
over a distance of the order of the range of /isr(l, 2) had already been noted in the pionneering 



(*)A different approach has been proposed recently by Alastuey and Ballenegger 
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work of Nienhuis and Deutch ^ . We stress that eq. (|12|) does not express the permittivity in 
terms of local fluctuations of the dipole moment, because hsr{l, 2) is only a part of the total 
correlation function. Whether such a relation can be established in the general case of a local 
permittivity tensor e(r) is under current study. 

Let us apply the result H12() to the situation illustrated in Fig. 1. The symmetry of the 
planar interface implies that e(r) = e{z) is diagonal, with exx — Syy — e" and ezz = e^- When 
the external field Eq — E^^e^ is parallel to the interface, V A E{r) — implies ^.©"(z) = 0, 
and hence E{z) = E^^ix is constant throughout space. In this particular case, eq. Hll|l follows 
exactly from eq. 0. Very close to the interface, the permittivity e"(2;) is expected to vary on 
a microscopic scale, and the concept of a local permittivity becomes really meaningful only 
for sufficiently large z, where e^^{z) varies over mesoscopic distances. 

In the case of a perpendicular external field E^^t = E-^Bz, the Maxwell field E{z) displays 
a non trivial profile along the z-axis. The local permittivity e-^{z) is therefore well defined 
only in regions sufficiently far from the interface, where the variations of E{z) are neglegible 
over the range of decay of the correlation function hsr{l,2). 

The leading deviations of the local permittivity (z) from the bulk dielectric constant e 
of the fiuid can be calculated explicitely for large z and low densities. For a planar interface, 
the interaction U{1) [eq. Q] of a dipole pi with its own image is given by 



C/(l) 



p2 1 - e' 1 + cos2 01 



2 1 



(2zi)3 



Pi ■ ez. 



(14) 



We consider first the case e' ^ I. The one-body density is then given at low density by 
p{l) ~ /ci/(47r)[l — PU{1)] when zi — > oo. Substituting in eq. p2|l and neglecting all terms 
of order p^, we see that the presence of the interface adds to Debye's well known low-density 
result e ~ 1 -(- the following surface corrections 



(^) 



l-e' 9y I3p^ 



1 - e' 33y 
1 + e' 80 z3 



(15) 



When e' = 1, U{1) vanishes and there is no surface term at order p. The correction of order 
p^ can be obtained exactly by expanding p(l) at low density. Treating the wall as a giant 
impenetrable particle (labelled 'w'), we have p(l) = p^'^\l,w) / p{w). The two-body density 
p'^''(l, w) can then be expanded using Mayer graphs. At order p^, a single graph contributes 
to p(l) when zi oo, and we find (see for the details of this calculation) 



.(1) . 



1 



16zf 



Oip') 



Zl 



oo. 



(16) 



Notice that the factor 3y/2 can be interpreted as the image charge factor (e — l)/(e -f 1), 
expected from a macroscopic electrostatics calculation, when e 1 + 3y. Inserting eq. 1)16(1 
into H12() shows that the deviations of €ij{z) to the bulk e are given in the case e' = 1 (up to 
order p"^ included) by the same expressions lfTK|l with additional factors 3y/2. 

The above low-density results agree with the asymptotic expression of e(z) obtained in [31 
using a generalization of the Onsager cavity model to the case of a planar interface. An 
exact asymptotic expression for p{l) at large distances for finite densities has been obtained 
by Badiali TU' and Zhang, Badiali and Su m] using Mayer diagrammatic techniques. Their 
general expression for p(l) involves a term with an integral over the three-body direct corre- 
lation function c(l, 2, 3). A similar analysis applied to Cy (z) shows that this function decays 
asymptotically as z^'^ , with a coefficient that involves integrals over two, three and four-body 
direct correlation functions. 
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Potential of mean force. - The screening properties of the polar fluid close to the interface 
can be probed by introducing two fixed ions of charge qi and (72 at ri = and r2 = 

{R2tZ2), with zi^Z2 > 0. Our aim is to determine the exact long-range behaviour of the 
potential of mean force w{\,'2,) = w{qi, q2, zi, Z2, R12) between the ions when their relative 
distance i?i2 along the wall is very large. From the Coulomb interaction energy between 
the two ions is 

Kc(l,2) . ,r,2G(n,r2) -^ + \-^^- (17) 

ki2| l + e' |t-i*2 

The interaction energy Vcd{i,j) = Vdc{j,i) between an ion i — 1,2 and a dipole j has a 
short-range part j) and a long-range part 

Vcdihj) =vcdii,j) + -—v*^{i,j), Wcd(i, j) = -(Pj ■ Vj)— ^. (18) 
1 + e \rtj 

The image interaction v*^{i,j) is obtained from Wcd(*,i) by replacing by r*y Similarly, the 
interaction energy Vddihj) between two solvent molecules has also a short-range part v^^{ijj) 
and a long-range part 

1 - e' 

Vddihj) = vddiij) + Y:j7^"dd(*>i) (19) 

where ?;dd(l,2) is defined in (|7|) and v*^^{i, j) is the image dipole interaction. 

The potential of mean force — /3w(l,2) is given formally by — /3T4c(l;2) plus the sum of 
all connected Mayer graphs made of two root points representing the ions, internal points 
of weight p{i) — p{zi,Pi) representing solvent molecules, and Mayer bonds / which take 
the value fddihj) = exp[— /3Vdd(*, i)] — 1 between two solvent molecules and fcdihj) = 
exp[—f3Vcd{hj)] — 1 between an ion i and a dipole j (see for example Q2|)- The graphs 
must be free of articulation points, and the root points must not form an articulation pair. 

The charge-dipole bonds fed and fdc decay as l/rjj at large distances. We decompose 
therefore fed (resp. fdc) into /cd(«, j) = /c'd(i, j) + </'cd(i, j), where (t>edihj) = -/3Wcd(«> j) is the 
electrostatic charge-dipole interaction and /^^ decays more rapidly. The graphs contributing 
to w{l,2) at large relative distances R12 are those which have only a single bond 4']^^{l,i) 
attached to the ion 1, and a single bond (p^dcij' 2) attached to the ion 2: 



-/3u;(l,2)^~/3Ke(l,2) + Q,'- + Q,'- \, i?i2^oo. (20) 

12 1 2 

In the second graph, the solvent molecules 3 and 4 are connected by the total (inhomogeneous) 
dipole-dipole correlation function /i(3,4). The contribution of the first graph is 



drg / dr!3 cj,':^{l, 3)p(3)^:fe(3, 2) = <^cd + -— ^d * + —— 03c , (21) 
23>o J 1 + e J L 1 -|- e J 

where we denote the convolution by a star. If the convolution 0cd * 4>dc were perfomed over 
the whole space {i.e. with /9(3) replaced by a constant), it is clear from dimensional analysis 
that it would lead to a term decaying as 1/r (because we would integrate two bonds decaying 
as 1 /r^ over a 3-dimensional space) . The regions where the point 3 is close to the root points 
leads to terms decaying as since (f)cd and (f)dc ^ . The asymptotic behaviour of the 

convolution is hence determined by the integration over the regions where the point 3 is far 
away from the points 1 and 2. 
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When focusing on the asymptotic behaviour R12 ^ 00 of the convolution (|21|l . it is legiti- 
mate to replace p(3) by the step profile 9(z3)p/(47r). Indeed, let 0cd*0dc = qiq2P^ J{zi, Z2, R12) 
and take the limit R12 ^ cxo by setting R12 — R12/1 with 7^0. From Ijf 8() . we have 



J{zi,Z2, ) 

7 



dr^ 



3 P3 



V' 



f 



P3 



V' 



1 



z'>0 



3 -^12 
1-72=12 
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zi,P3) (22) 



where we introduced the change of variable r'^ — j{r^ — ri). J{zi, Z2, R12/"/) is hence pro- 
portional to 7 as 7 ^ (assuming that Z12 does not increase faster than I/7), so that 
J{zi,Z2, R12) decays as I/R12 at large distances. Moreover, hm^^o Piz's/j + zi^p^) = p/(47r), 
and the convolutions in eq. I^T)l can therefore be calculated with the step profile 9(z3)yo/(47r) 
in place of p{3) when R12 — > 00. The results for these integrals are shown on Table I. 

Table I - Convolutions over the region z > of interaction bonds between dipoles and charges: 



/ 




0dd 


(pdc 


<t>dc 




-3j/0cd + ^0cd 






2 V^cc 


<^cd 


-f0c*d 


-'f^cd 


2 r^cc 





When the decomposition c(l,2) ~ Csr(l,2) + 0dd(l,2) + (p'^^{l,2) is introduced in the 
Ornstein-Zernike equation, the correlation function h{l, 2) splits into the short-range part 
hsi(l, 2) [see eq. ()10|) ] and a long-range part h\^{l, 2). Consider the second graph of H20(l with 

in place of h and add this term to eq. 



Z3>n 



d3d4</)|,'d(l,3)£[(5(3,4) + £-/^,,(3,4)j0[;,(4,2), i?i2 ^ 00. (23) 



Since the points 3 and 4 are close together (because /isr(3,4) is short-ranged) and far away 
from the root points, they can be treated as a hypervertex. Using 1)18(1 and p3|l . expr. H23|l 
becomes 



e - 1 



d3^^'d(l,3)^</.^'e(3,2)- 



e - 1 r 



Z3>0 



4tt 



0cd 



l-e' 
1 + e' 



(t>dc 



l-e' 
1 + e' 



(24) 



The long-range part /iir(3,4) has been studied by Badiali ^3] and Badiali and Forstmann |14j . 
For points 3 and 4 far apart and far away from the interface, it is given by 



/ilr(3,4) 



e - 1 \ 1 
e 



32/ 



0dd(3,4) 



:7<^Sd(3,4) 



The second graph in H2U|) with /iir in place of h gives thus the contribution 



3j/ 



l-e' 



-1 V^cd 



0dd(3,4) + ^0*,(3,4) 



fdc 



(25) 



(26) 



According to l|2U|) . the potential of mean force is given at large distances by V^c minus the 
expressions (|24f) and H2()l) over /3. Using the results shown on Table I for the convolutions, we 
find after some simple algebra: 



V. BALLENEGGER et al: LOCAL DIELECTRIC PERMITTIVITY NEAR AN INTERFACE 



7 



The present calculation shows therefore that, even though the ions are close to the interface, 
their effective interaction at large relative distances along the wall is given by the standard 
formula of macroscopic electrostatics. Local effective dielectric functions for solvated macro- 
molecules obtained from continuum electrostatics models have thus some microscopic 
basis. 

* * * 
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